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ABSTRACT 

The globular cluster Gl (May all II) in M31 is the most massive (~10 7 M Q ) stellar cluster in the Local 
Group, and it has the highest central velocity dispersion (^28 kms~ ). It has been claimed to host a central 
~20,000 Mq black hole, but these claims have been controversial. We report here the XMM-Newton detection 
of X-ray emission from Gl at a level of L x w 2 x 10 36 erg s" 1 . This emission could be the result of Bondi-Hoyle 
accretion of ionized cluster gas by a central black hole, or it could be produced by a conventional low-mass 
X-ray binary. A precise localization of the X-ray emission, which is not possible with the current XMM data, 
could distinguish between these possibilities. While such a measurement may be difficult, it is of sufficient 
potential importance to pursue. 

Subject headings: globular clusters: individual (Mayall II = Gl) — X-rays: binaries 



1. INTRODUCTION 

Dynamical modeling of the kinematic data of numerous 
nearby galaxies has revealed that central dar k objects of 10 6 - 
10 10 Mf?> are nearly u biquitous (e.g- lKormendv & Richstonel 
1 1 9951 iMagorria net alJI 1 9981) . These central dark objects are 
generally assumed to be supermassive black holes, the rem- 
nants of quasars which were much more numerous in the past. 
Another population of black holes, those with masses of ~3- 
20 M Q , is known from determinations of the binary param- 
eters of 18 X -ray binaries in the Milky Way and Magellanic 
Clouds (e.g., Remil lard & McClintockl 120061) . Bridging the 
gap between these two populations, i.e., finding intermediate 
mass black holes (IMBHs), has been a long-standing goal. 

Some evidence for IMBHs has come from application 
of the dynamical modeling technique to globular clusters, 
which has found ce ntral dar k objects in a couple of cases. 
iGerssen et alJ (12002) and lvan der Marel et alJ J200 2) claimed 
that the Milky Way globular cluster M15 hosted a cen- 
tral IMBH with a mass of (3.9 ± 2.2) x 10 3 Mq, and 
iGebhardt et alJ (12002b argued that the Andromeda cluster Gl 
hosted an IMBH of 2^1 x 10 4 M W . Both of these claims were 
challenged by Baumgardt et al. (2003afl). 

The case for M15 hosting a central IM BH was we akened by 
the discovery of a mislabeling of axes (Dull et al. 2003), but 
the case for Gl has recently been s trengthened by new, high- 
quality data (Gebha rdt et al.l l2005). So me relevant physical 
parameters for Gl are given in Tabled Mevlan et al. (2001) 
have suggested that Gl could be the stripped core of a dwarf 
elliptical galaxy. 

If a globular cluster does indeed host a central IMBH, and if 
it contains some amount of ionized gas (as has been observed 
in 47 Tuc: lFreire et al.l2001l) . one might reasonably expect the 
IMBH to accrete the gas via a Bondi-Hoyle-Lyttleton process 
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llHovle & Lvtfletonll941tlBondi & Hovlell944tlBondil 19521) 
and become an observable X-ray and/or radio source (e. g . , 
lGrindlavll978tlPfahl & Rappaporl200U lMaccaronel2 004 ) . 

iHoetaLl J2003I) reported that a search for a central X-ray 
source in a Chandra observation of M15 yielded a null result. 
We report here on the detection of X-rays from G 1 with XMM- 
Newton. This X-ray detection may be emission from a central 
IMBH but is also consistent with emission from a low-mass 
X-ray binary (LMXB). The observations are described in §|2 
and the analysis in §[3] We discuss the possible interpretations 
in §|4]and briefly conclude in §|5] 

2. OBSERVATIONS 

XMM has observed Gl on three occasions. A log of the 
observations is given in Table |3 We downloaded the re- 
processed European Photon Imaging Camera (EPIC) data for 
each observation from the XMM Science Archive 2 . EPIC con- 
sists of three separate cameras: one PN CCD array and two 
MOS CCD arrays. The first observation had Gl at the aim- 
point, but the other two included Gl at a considerable off-axis 
angle, which reduced the effective area due to vignetting. To 
give a sense of this, we list in Table |2 the effective area of 
the PN instrument at 1 keV, as determined from the ancillary 
response files. The XMM Optical Monitor was off for the first 
observation and did not include Gl in the field of view for the 
other two observations. 

We used the XMM Science Analysis System (SAS) version 
6.5.0 (obtained from the XMM Science Operations Centre 3 ) 
for the reduction of the data. The second observation suffered 
from intense background flaring, and we filtered out the worst 
flaring periods to use only a few timeframes which were rel- 
atively unaffected. This resulted in a significantly reduced 
exposure time. 

We extracted spectra from the location of Gl for each in- 

2 |http : / /xmm . vilspa ■ esa ■ es/external/xmm_data_acc/xsa| 
3 |h^T^^7^gnn^^sa^^esa^^^7r 
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FIG. 1 . — Left: An 8' X 8' XMM-Newton image of Gl with photons in the range 0.5-1.2 keV shown in red, those in 1.2-2.5 keV shown in green, and those in 
2.5-6 keV shown in blue. The image has been smoothed by a Gaussian with a full width at half maximum of 12" and is displayed with a high contrast ratio in 
order to suppress the background. The small box indicates the area of the HST image shown on the right. Right: A 20" X 20" HST image of Gl taken with the 
F555W filter. 



strument for each observation and constructed response ma- 
trices and ancillary response files using SAS. We used large, 
nearby, source-free regions to extract background spectra. 
The photometry based on these extractions is given in Table[3] 
We have combined the detected counts from the PN and both 
MOS instruments. Because the exposure time and effective 
area vary greatly between instruments and observations, we 
report both net counts and net count flux. We use the mean 
effective areas over 0.2 to 10 keV to calculate the fluxes. 

An image of the first observation is shown in Figure [2 in 
which the photons have been color-coded according to energy. 
We also show a Hubble Space Telescope image of Gl taken 
with the Advanced Camera for Surveys with the F555 W filter. 
This multi drizzled image was downloaded from the Multimis- 
sion Archive at the Space Telescope Science Institute 4 . 

3. ANALYSIS 

Since the count fluxes of all three observations are consis- 
tent with each other, we use only the first for spectral analysis 
because it has the best signal to noise ratio. We jointly fit the 
three spectra (from the PN and both MOS instruments) from 
the first observation in Sherpa ( Freeman et al.l200ll) using an 
absorbed power law model. We considered only the 0.2-10 
ke V band, and the spectra were grouped to c ontain at least 3 
counts per bin. We used \ 2 statistics with the Gehrels ( 1986) 
estimate of the standard deviation of each bin. The column 
density was constrained to be at least equal to the G alactic 
value of 6.57 x 10 2() cm -2 (Dickey & Lockman 1990). The 
parameters of the best fit are given in Table El They are not 
well constrained. The unabsorbed 0.2-10 keV luminosity is 
~2xl0 36 ergs _I . 

4. DISCUSSION 

4 |http : / /archive . stsci . edu/| 



4.1. IMBH accreting cluster gas 

A black hole of mass Mbh, moving with relative speed v 
through a gas of ambient density p and sound speed c s , nom- 
inally accretes at the Bondi-Hoyle-Lyttleton rate (Hoyle & 
Lyttleton 1941; Bondi & Hoyle 1944; Bondi 1952; hereafter 
BHL) of: 

M~4^(GM BH ) 2 p(v 2 + C 2 r 3 / 2 . (1) 

If we define V = (v 2 + c 2 ) I//2 and p = m p n, where m p is the 
proton mass and n is the hydrogen number density, then for a 
black-hole mass of 10 4 M Q the corresponding X-ray luminos- 
ity is given by: 

L x ~ eMc 2 ~ (2) 
SxlO 37 ^) 2 ^)^^)^ (3) 

where we have taken e = 0.057 for the efficiency of converting 
mass to radiant energy in the inner disk of a Schwarzschild 
BH, and V has been normalized to that expected for a gas 
at T = 10 4 K with no significant bulk motion relative to the 
BH. (Of course, the fraction of the energy emitted in the X- 
ray band may be considerably smaller than unity.) In order 
for hydrodynamic accretion to take place, the gas most likely 
has to be ionized (for mean free paths to be shorter than the 
accretion radius), and an ionized state may be achieved by 
the UV from only one or two post-AGB stars in the cluster. 
IPfahT fe Rappaport (2001) argued that the wind from cluster 
AGB stars could maintain the intercluster gas density at ~0. 1— 
1 cm" 3 , even after the cluster is completely stripped of its gas 
upon passing through the galactic plane of the host galaxy. 

Direct empirical evidence for a substantial ionized gas den- 
sity in 47 Tuc, via observations of differential dispersion mea- 
sures of the resident radio pulsars (iFreire et alJ2001l) . yields a 
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TABLE 1 
Properties of G 1 in M3 1 


Parameter 


Value 


Ref. 


Total Mass (M) 


1.5 X 1O 7 M 


[1] 


Central Density (po) 


4.7 X lO 5 M pc- 3 


[1] 


Central Velocity Dispersion (tro) 


27.8 km s" 1 


[1] 


Core Radius (r c ) 


0.14" 


[1] 


Half-Mass Radius (r h ) 


3.7" 


[1] 


Tidal Radius (r t ) 


54" 


[1] 


Distance (d) 


750 kpc 


[2] 



TABLE 2 
XMM-Newton OBSERVATIONS OF Gl 



REFERENCES. — m : MevlanetalJ<200ll) : [2] Freedm an et alj<200l) 



value of n = 0.067 ± 0.015 cm" 3 which confirms the possibil- 
ity that globular clusters may indeed contain a non-negligible 
density of ionized gas. 

The case of Bondi-Hoyle accretion by an IMBH is com- 
plicated by the large number of field stars that are likely to 
be contained within an accretion radius. The accretion radius 
used to derive eq. (1) is: 



2GM BH 



0.4 



V 



pc (4) 



V 2 "" \10 4 M Q J Vl5km/s, 

According to Fig. 6 o f lGebhardtet ail (12005). there are more 
than 10 stars within ~0.4 pc of the center of Gl. Thus, we 
should consider both the dynamical effect of such stars on 
any gas present and the direct injection of gas into the accre- 
tion volume via the stellar winds of these stars. We start by 
considering the collective effects of the stellar wind from the 
stars neighboring the IMBH. Each star has a stellar wind in- 
tensity and speed that depend sensitively on its evolutionary 
state. We assume that these winds collide and their bulk ki- 
netic energies are thermalized. 

In order to estimate the net mass loss and mean thermal ve- 
locities of the colliding winds, we consider the evolution of a 
single star near the turnoff mass (i.e., with M ~ 0.85M Q ). We 
take the evolution of this one star to represent an ensemble 
of stars in different evolutionary state s. We ut ilized a sim- 
plified evolution code (Soker & Rappaport 2000) with a wind 
loss prescription that is a modified version of Reimers' (1975) 
expression. The Reimers wind is enhanced as a function of 
the star's mass and radius (see eq. [15] in Soker & Rappaport 
2000) so as to reproduce the initial-final mass relation for the 
production of white dwarfs in single stars ( Weidemann fl993l 
We take the terminal wind speed, v w , to be equal to the escape 
speed from the stellar surface. During the evolution of our 
representative cluster star we tabulated two integrals: 



AKE = 



AGB 



1 ■ , 

-Mvldt ; AM = 



ZAMS 



AGB 



Mdt 



(5) 



ZAMS 



From these we compute the mass-weighted mean specific ki- 
netic energy and the mean mass-loss rate from the winds: 



< K, > 



<M> 



AKE /AM 
AM/ At , 



(6) 
(7) 



where At is the evolution time from the ZAMS to the end of 
the AGB phase. 

From these calculations we find the rms wind speed after 
collisional thermalization to be v w .rms = v t h ~ 55 km s" 1 , while 
the average value of M ~ 2 x 10~ u Mq yr" 1 . Of course, the 



Date 


ObsID 


Exp. Time (s) 


PN eff. area 






PN 


MOS 


(cm 2 at 1 keV) 


2001 Jan 11 


0065770101 


4486.5 


7325.1 


894.2 


2002 Dec 29 


0151580101 


1954.4+ 


3816.9 1 " 


482.1 


2003 Feb 06 


0151581101 


7271.3 


9563.6 


556.4 



+ Filtered for background flares. The total exposure times were 6106.8 s for 
the PN and 8189.4 s for both MOS cameras. 



bulk of the wind loss contribution comes from stars on the 
giant branch or AGB. If we utilize this wind speed in eq. (4) 
we estimate /? acc ~ 0.06 pc. Based on the central density po, 
there are ^425 Mq of stars in this volume, the bulk of which 
should be composed of the most massive cluster stars, near 
the turnoff or in the giant phase, which will have preferen- 
tially sunk to the cluster center. These ~500 stars could pro- 
duce a total mass injection rate, via stellar winds, >10~ 8 M Q 
yr" 1 . This is approximately 100 times more than is needed to 
power the observed L x from the vicinity of Gl. Finally, we 
point out that since the rms speed of the cluster stars is lower 
than the rms thermal velocities that we estimate (v t h ~ 55 km 
s" 1 ), the dynamical effects of the stars on the gas are probably 
negligible. 

4.2. LMXB(s) in outburst 

In the Milky Way, we know of 13 LMXBs in twelve 
globul ar clusters with L x > 10 36 ergs" 1 (for a review, see 
IVerbunt & Lewinl20b 6. and references therein). Five of these 
sources are recurre nt transients, and the other seven are per- 
sistently luminous. Pool ev et alJ 112003) have shown that the 
number of globular cluster close binaries observable as faint 
X-ray sources (e.g., LMXBs in quiescence, cataclysmic vari- 
ables, active main-sequence binaries) scales as the "encounter 
frequency" of the cluster, which goes roughly as po 15 r c 2 . In 
their units, Gl has an encounter frequency of ^7500 and 
would be expected to host ~75 LMXBs, the vast majority in 
quiescence. If some small number of those LMXBs were in 
outburst during each of the three XMM observations, it could 
explain the observed X-ray emission. 

4.3. Distinguishing the two possibilities 

Because the observed value of L x is compatible with ei- 
ther scenario, we discuss two possible ways to distinguish 
a central IMBH from an LMXB. The first is by its loca- 
tion with respect to the cluster center, where the more mas- 
sive the object the closer it is likely to be found to th e cen- 
ter (e.g., Jernig an & Clarldll979t lG"rindlav et al. 1984). Near 
the cluster core, the stellar density approaches a constant 
value, and the cluster potential grows as r 2 . If black holes 
(BHs), neutron stars (NSs), and field stars are in rough en- 
ergy equipartition, then we expect the rms black hole ve- 
locity to scale as vnsv^ns/^bh where M N s and v N s are 
typical neutron star masses and velocities. Thus, we would 
anticipate that any IMBH would be confined within a ra- 
dius of Jns v /Mimbh — 0.02 r^s, where r^s is the ra- 
dius containing about half the neutro n stars in the cluster (for 
a more detailed discussion see, e.g., Ijerni gan & Clarill97§ 
iBahcall & Wolil976l) . According to lHeinke et alJ J2003I) . the 
neutron star distribution in clusters follows that of the field 
stars, but with a characteristic radial scale that is ^60% that 
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TABLE 3 

X-ray Photometry of Gl in 0.2-10 keV 



Obs. Net Counts 


Net Flux (cnts cm 2 s 1 ) 


2001 Jan 11 51.5 ±8.6 

2002 Dec 29 21.9 ±14.4 

2003 Feb 06 20.0 ± 12.8 


(1.02±0.17) X 10~ 5 
(1.51 ±0.99) X 10~ 5 
(0.48 ± 0.31) x 10 -5 


TABLE 4 

Absorbed Power Law Fit Parameters 


Parameter 


Value 


Column density (lin) 
Photon index (a) 
Normalization at 1 keV 
X 2 /d.o.f 
Absorbed F x t 
Unabsorbed Fx' 
Unabsorbed L x ' 


6.6 +82 x 10 20 cm- 2 

1 '^-0.6 

2.4*1 x 10 ^ cntskeV" 1 cm" 2 s" 1 
4.5/16 

2.8 X 10~ 14 erg cm -2 s - ' 
3.0x 10~ 14 ergcnrV 
2.0 xlO 36 ergs" 1 


+0.2-10 keV 



of the field stars. Therefore, we expect half of the LMXBs 
in Gl to be within ~0.6r/, ~ 2'.'2 of the center, but a 20,000 
Mq IMBH to be within 50 milliarcsec of the center. Unfortu- 
nately, the spatial resolution of XMM does not allow for such a 
test, but a Chandra observation could easily localize the X-ray 
emission to a fraction of an arcsecond. If there are at least a 
few X-ray /optical coincidences in the Chandra field of view, it 
should be possible to register the Chandra and HST frames to 
0"l-0"2 accuracy. With such a registration, an IMBH should 
be essentially coincident with the center of Gl, and there is a 
good chance that any detected LMXB would be certifiably re- 
moved from the center. The probabilities that a cluster LMXB 
would lie within 0."3, 0"2, and 0"1 of the center are 0.16, 
0.11, and 0.04, respectively. 

Although accretion disk spectra are notoriously difficult to 
calculate from first principles, an IMBH and an LMXB may 
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have observably different spectra. It has been suggested (e.g., 
iKaaret et ail 120031 iMiller et al1l2003l) that a cool multicolor 
disk spectral component, indicating a relatively low inner disk 
temperature, might indicate the presence of an IMBH. The 
sense of the effect is that the larger the mass of the BH accre- 
tor, the lower the temperature, T, of the inner edge of the disk, 
which scales as: 

T t cx (M B hM/V, 3 ) 1/4 K M i/4 Mbh /2 (8) 

for a simple thin-disk model. Since some Galactic BH 
sources, with M Bll ~ 8M Q and M ~ 10" 8 M© yr" 1 have 7] ~ 
1 keV, then an IMBH in Gl might have T, ■ ~ 10 eV. Thus, an 
association of the X-ray source with an IMBH, as opposed to 
an ordinary LMXB containing a neutron star or stellar-mass 
BH, could be made on the basis of a very soft observed spec- 
tral component. We know of no reported observations of such 
a component in Gl . 

5. CONCLUSIONS 

We have reported on the detection of X-ray emission from 
the massive star cluster Gl in Andromeda, observed three 
times with XMM. There are two very different possibilities 
for the origin of the X-rays: an IMBH accreting cluster gas 
or one or more LMXBs in outburst. Unfortunately, the XMM 
data do not allow for a discrimination of these two scenarios. 
Future observations may reveal that the source of the X-rays is 
displaced from the center of the cluster (indicating an LMXB) 
or found at its center (indicating an IMBH), and they may de- 
tect a very soft spectral component (supporting the presence 
of an IMBH). 
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